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AnSTRACT 

Diffusion  coefficients  of  four  binary  electrolyte 
systems  vjere  measured  using  the  capillary-reservoir 
method.  The  four  systems  are:  1.0  molar  KCl  in  aqueous 
solution  diffusing  into  H2O  ; 0.01  molar  ZnO  in  40/o  KOH 
solution  diffusing  into  kO^  KOH  ; O.05  molar  aqueous  CuSO^^ 
diffusing  into  H^O  ; and  a series  of  0.005  - 0.10  molar 
CuSOji^^  in  0.5  molar  H^SO^  diffusing  into  0.5M  H£S0f.  The 
ZnO/KOH  and  CuSO^/H^SO^  systems  are  of  current  interest 
to  electrochemical  engineers,  and  demonstrate  the  use  of 
this  technique  for  binary  electrolytes. 

The  diffusion  coefficient  of  the  ZnO/KOH^ystera 
at  25®C  v;as  measured  as  7.22  ± I.90  x 10“^  cm^/sec. 

Across  the  range  of  concentrations  of  CuSOji^  in  H2S0jij^ 
described  above,  the  average  diffusion  coefficient  v/as 
measured  as  5.^3  i 0,59  x 10-^  cm“/sec  at  25°C,  and 
6,675  - 0,97  X 10“^  cm^/sec  at  35.5®0,  A slight  effect 
of  diffusant  concentration  on  the  diffusion  coefficient 
was  found. 

The  CuS0|^/H20  system  was  studied  in  two 
different  versions  of  the  capillary  technique:  one 


I 


in  which  the  reservoir  v;as  stirred,  and  the  other  under 
no-stir  conditions.  Results  of  these  experiments 
Indicate  that  the  hydrodynamic  effect  caused  by  stirring 
leads  to  an  apparent  diffusion  coefficient  higher  than 
would  be  expected,  unless  sufficient  diffusion  tine 
is  allowed.  Present  data  lead  to  the  conclusion  that 
the  no-stir  model  is  the  most  practical  choice  for 
obtaining  the  diffusion  coefficient  of  this  and 
similar  binary  electrolyte  systems. 
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Electrochemical  systems,  such  as  electric  batteries 
and  electroplating  tanks,  cannot  be  v.’ell-designed  without 
knovrledge  of  the  transport  properties  of  the  system,  such 
as  the  diffusion  coefficient,  viscosity,  and  the  electrical 
conductivity.  For  example,  the  relation  which  gives  the 
limiting  current  density  at  an  electrode  in  a given  electrolyte 
is  a function  of  the  diffusion  coefficient  of  the  migrating 
ionic  species  in  supporting  electrolyte. 

For  many  binary  liquid  pairs,  fluid  properties  ouch 
as  density  and  viscosity  are  tabulated  as  functions  of 
composition,  temperature,  and  pressure.  Where  not 
tabulated,  these  properties  may  often  be  very  closely 
estimated  with  well-developed  approximations  found  in 
scientific  literature.  Unfortunately,  though,  the 
diffusion  coefficients  of  many  electrolytes  are  not 
found  in  the  literature.  It  is  the  intention  of 
this  study  to  explore  some  aspects  of  a method  of 
measuring  diffusion  coefficients  of  aqueous  electrolytes. 

If  this  method  can  be  shovm  to  be  reliable  and  accurate. 

It  will  be  possible  to  accurately  determine  the  diffusion 
coefficients  of  electrolyte  systems  not  yet  discussed 
in  the  literature. 


Most  analytical  techniques  for  measuring  a property 
of  a system  are  based  on  a theoretical  model  of  the 
system  under  consideration.  If  the  model  is  an  accurate 
representation  of  the  system,  then  the  value  of  the 
measured  property  is  believable.  For  example,  a 
thermometer,  which  is  thought  of  as  a device  for 
measuring  temperature,  represents  a model  of  the 
relationship  betvieen  the  expansion  of  mercury  and 
a change  in  temperature,'  By  calibrating  the 
thermometer  at  two  standard  (reference)  points, 
such  as  the  melting  and  boiling  points  of  v:ater 
at  one  atm,'  pressure,  it  becomes  a standardized, 
reproducible  measurement  device.' 

Similarly,'  the  techniques  used  for  measurement 
of  liquid-liquid  diffusion  in  a binary  system  arc 
based  on  a model  of  Pick's  First  Law  of  Diffusion.' 

The  capillary  method  of  diffusivity  measurom.ent 
represents  one  such  model,'  In  the  development  of 
the  model  used  in  the  capillary  method,  boundary 
conditions  must  be  chosen  for  the  solution  of  the 
diffusion  equation,'  In  order  to  obtain  a simple 
analytical  solution  of  the  equation,  it  is 
necessary  to  maintain  a constant  zero-concentration 
of  diffusant  at  the  capillary  mouth.  Unfortunately,' 


doing  GO  cauGGs  ujii^anted  side  effects  due  to  fluid  | 

flowing  past  the  capillary.’  And  if  concentration  of  '■ 
I diffusant  at  the  capillary  is  not  kept  at  zero,  the 

solution  of  the  diffusion  equation  becomes  very 
complicated  and  must  bo  approximated  by  numerical 
methods,’  In  this  study,  experiments  based  on  both 
models  were  done,  and  compared,  vzith  the  intention 
of  being  able  to  make  a recommendation  for  future 
experimentation  v^ith  the  capillary  method.’ 

The  electrolyte  systems  studied  in  this  series 
of  experiments  include;  l.'O  molar  potassium  chloride 
in  water  diffus^-ng  into  vjater  at  25®C.  ; O.'Ol  molar 
zinc  oxide  in  (wt.)  potassium  hydroxide  diffusing 
into  KOH  at  25®C,  ; 0,05  molar  cupric  sulfate 
in  water  diffusing  into  26®  water  ; and  a series 
of  lovr-concentration  (less  than  or  equal  to  0,1  molar) 
solutions  of  cupric  sulfate  in  half-molar  sulfuric 
acid,'  diffusing  into  0.5H  H2SO4  at  both  250C  and  35°C.' 


II.  Hinf-orlcaT.  Backnr'oimd 


A.  Hlntory  of  Dlffuslvitv  Koanureniont 
1.  Methods  of  Diffasivloy  ^!oasuron:^■:nt 

There  exist  many  techniques  for  measurement  of 
the  diffusivity  in  a binary  system,  all  of  vdiich  have 
the  same  basic  pattern:  Initially  species  "A"  is  in 
one  container,  isolated  from  another  container  holding 
species  “B",  For  a measured  time  period,  the  tvio 
species  contact  each  other,  then  the  two  are  separated. 

The  contents  of  each  container  afo  analysed  in  order 
to  determine  how  much  "A”  is  in  the  container  v.'hich 
originally  held  only  "3",  and  likewise,  how  much  "B" 
is  in  the  "A"  container.  The  textbook  Tec’inioues  of 
Orrranlc  Chemistry  (1)  contains  a chapter  called  "Diffusivity 
in  which  many  of  these  measurement  techniques  are 
discussed.  For  example,  it  is  sometimes  useful  to 
have  the  two  diffusion  cells  cylindrical  in  shape 
and  allov;  contact  through  one  open  end.  This  technique 
is  often  used  in  measuring  the  diffusivity  (the  diffusion 
‘ coefficient)  of  binary  gas  mixtures.  A variation  of 
this  method,  in  which  the  two  cylinders  arc  separated 
by  semi-permeable  membrane,  allov;s  measurement  of 
the  transport  of  gases  or  liquids  through  such  a 
membrane . 


, In  measuring  the  diffusion  coefficients  of  electrolytes, 

i it  is  useful  to  measure  the  conductance  of  the  electrolyte 

I 

i 

! before  and  after  diffusion.'  In  The  Physical  Chenistry 

■ o£  Electrolytes  (2)  , this  technique  is  developed  in 

, detail,  and  diffusion  coefficients  for  nany  aqueous 

I electrolytes  are  tabulated.  Experimental  data  based 

* 

1 dLa-to. 

on  Hamed’s  conductometric^ may  be  found  in  contemporary 
literature  describing  diffusivities  of  very  dilute 
aqueous  salt  solutions  (3)  . 

Other  established  diffusion  measurement  techniques 
•used  with  aqueous  electrolytes  include  measuring  the 

weight  change  of  a diffusant-f illed  oorous  disc  (4)  . 

i 

i and  measuring  the  conductance  of  a diffusing  system 

with  a rotating  disc  electrode  (5)  .' 

Z,'  The  Capillary-Reservoir  Method 

Introduced  in  the  late  19^0's,  the  capillary 
method  of  diffusivity  measurement  remains  one  of  the 
most  versatile  and  reliable  techniques  inuse.'  As 
originally  described  by  Anderson  and  Saddinton  (6)  , 
a capillary  tube  with  small  diametcr-to-length  ratio 
is  filled  with  a kno’-m  am.ount  of  diffusant  ("A")  and 
is  placed  into  a container  (reservoir)  filled  with 
"B",  for  a measured  time  period.  After  the  diffusion 
time  has  elapsed,  the  capillary,  is  removed  and  its 


contents  are  analyzed  in  order  to  determine  hov;  much  "A" 
remains,  Anderson  and  Saddington  used  radiotracers 
to  study  the  self-diffusion  of  some  ionic  species  in 
solution,' 

Many  applications  of  this  technique  have  been 
reported,'  Radiotracers  v/ere  used  in  studying  diffusion 
in  liquid  metals  (7)»(8)  , and  with  molten  salts  (9),  (10), 
The  radiotracer  method  of  analysis  has  been  used  v;ith  the 
capillary  method  in  studying  the  diffusion  coefficients 
of  some  organic  liquids  (11)  , and  of  liquified  inert 
gases  (12)  ,'  Because  of  the  intrinsic  simplicity  of 
the  capillary-reservoir  technique,  it  could  conceivably 
be  used  to  e'valuate  diffusion  coefficients  in  any  binary 
liquid  system,  as  long  as  some  method  exists  vrith 
which  to  determine  the  contents  of  the  capillary  before 
and  after  the  allotted  diffusion  time  has  elapsed. 

There  are  nany  features  of  the  capillary-reservoir 
technique  which  make  it  desirable  for  use  in  measuring 
diffusivities  of  a large  variety  of  liquid  pairs. 

The  apparatus  need  not  be  large;  the  size  of  the  capillary 
is  limited  only  by  the  amount  of  diffusant  necessary 
for  accurate  determination  of  capillary  contents. 

Because  the  system  may  be  kept  quite  small,  it  is  feasible 
to  conduct  diffusion  experiments  at  extreme  temperatures 


and  pressures  (13)»(l^) 


It  is  advantageous  to  minimize  capillary  length  for 
another  reason:  the  experimentally-determined  diffusion 
coefficient  is  a function  of  the  ratio  of  f inal-to-initial 
concentration  of  diffusant  in  the  capillary.  Thus,  for 
I a given  diffusion  time  period,  greater  accuracy  is 

possible  with  shorter  capillaries.  (See  Appendix  "B"). 

Despite  these  positive  features,  the  capillary 
technique  is  not  used  commonly  because  of  disagreement 
among  users  of  the  technique  v;ith  regard  to  the 
choice  of  boundary  conditions  for  the  diffusion  equation. 
There  are  tvfo  mathematical  models  being  discussed  by 
researchers:  one  is  called  the  "stirring"  model,  and 
the  other,  "no-stirring". 

lai_  The  Stirring:  Model 

During  diffusion,  the  capillary  sits  in  a liquid 
filled  reservoir.  In  order  to  achieve  a complete  analytical 
solution  of  the  diffusion  equation,  it  is  necessary  to 
maintain  the  capillary-reservoir  interface  (the  mouth  of 
the  capillary)  at  zero-concentration  of  diffusant.  This 
has  been  accomplished  in  the  past  by  stirring  (flowing) 
the  reservoir  fluid  rapidly  enough  to  remove  diffusant 


from  the  interface  as  soon  as  it  leaves  the  capillary. 
Unfortunately,  this  stirring  causes  a "sweeping  out" 
of  diffusant  from  belovr  the  capillary  mouth  (15). 

The  region  in  v/hich  this  "sweeping  out"  occurs  has 
been  termed  the  "Aa  -region",  since  mass  transport 
due  to  diffusion  only  occurs  through  a reduced  length 
of  capillary.  The  existence  of  this  region  has  been 
confirmed  in  the  laboratory  (16)  , and  the  magnitude 
ofA^  has  been  measured  by  photographic  methods  (15). 

A recent  paper  (17)  describes  computer  modeling 
of  the  streamlines  associated  with  flov7  past  the 
capillary  mouth,  so  as  to  allov;  estimation  of  the 
magnitude  of  A . This  model  is  based  on  work 
describing  flow  of  a viscous  liquid  past  a sharp 
corner  (18).  A pictorial  representation  of  this 
phenomenon  is  shovm  in  Figure  la.  In  Figure  lb,  the 
effect  of  viscous  forces  acting  on  the  liquid  is 
seen,  as  the  streamlines  become  horizontal,  leading 
out  of  the  capillary.  Figure  Ic  shows  the  overall 
A ^-effect  near  the  capillary  mouth.  Using  potassium 
permanganate  (deep  purple  in  color)  as  the  diffusant, 
and  water  in  the  reservoir,  Nanis,  Bockris,  and 
Richards  photogi'aphcd  the  A ^-effect  in  a stirred 
capillary-reservoir  system.  One  of  their  photographs. 


Pifcure  1 . Schematic  representation 


j of  the  flow  patterns  in  the  capillary 

' •.  leading  to  the  development  of  the 

' _ ■ -region,  (18). 


shovrins  diffuoant  streaming  out  of  the  capillary 
from  below  the  capillary  mouth,  is  shovm  in  Figure  2, 

Nanis  e:t  (15)  give  the  magnitude  of  the 
-region  as  a function  of  stirring  rate.  Calculation 
of  a Reynold's  Number  (Wro)  was  based  on  the  rotation 
rate  of  the  capillary  (cu) , the  radius  of  the  rotating 
arm  (r),  the  kinematic  viscosity  (v),  and  the  inside 
diameter  of  the  capillary  (d) , leading  to  the 
following  expression: 

Nge  * 2 7TrWd  /V  (2.1) 

Against  the  Reynold's  Number  they  plotted  a dimensionless 

liial  is,  ojC  divided  by  tne  inrernai  aiamecer  or 
the  capillary  (d).  This  plot,  shcvni  in  Figure  3, 
indicates  that  increases  v?ith  stirring  rate, 

approaching  a value  of  unity  at  large  The 

solid  curve,  of  Figure  3 is  an  approximation  which 
fits  the  data  at  lov;  Reynold's  Numbers  (less  than  10). 
This  approximation  may  be  expressed  as: 

A^/d  = 0.17Ngg  - O.Cl3Ngg2  + 0.00055Ngg^  (2.2) 

Experiments  at  high  Reynold's  Numbers  (Re  i 100) 
were  done  v;ith  liquid  metals  (19)  and,  when  combined 
with  the  data  in  Figure  3»  suggest  a value  of  A^/d 
approaching  zero  at  low  Reynold's  Numbers.  Because  of 


Photograph  of  KKrO^,  diffusing  into  water 
in  a stirred-reservoir  system  (15).  The 
effect  of  the  stirring  can  he  easily  seen 
in  the  capillary  on  the  right  (in  focus),' 


Plfrure  3.  Plot  of  Reynolds  Number 
against  A5/d  for  various  sizes  of 
» capillaries,  (15) 
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the  small  number  of  data  points  for  very  high  Reynold's 
Numbers,  no  significant  statement  can  be  made  about  the 
-effect  in  that  region.  In  that  study,  hov;eyer,  a 
relationship  was  suggested  betvfeen  A,J2  and  that 
could  be  used  as  a guide  in  estimating  the  expected 
error  in  diffusivity  calculations  due  to  the  -effect: 

= 12  (I®-)  (2.-3) 

Previous  experimentation  v/ith  the  A<^-effect  (20) 
has  implied  that  the  magnitude  of  the  A /-region  m.ay  be 
time-dependent.  Experiments  v.’ere  done  for  long  tim.es 
(several  days)  and  for  short  tir.es  (a  few  hours),  all 

the  apparent  diffusivity  for  short-time  experiments 
Vfas  consistently  higher  than  the  apparent  diffusivity 
for  long-time  experiments.  This  is  probably  due  to 
the  effect  of  the  ^region  growing  v:ith  time  until 
it  reaches  a constant  value,  after  which  the  diffusion 
effect  becomes  dominant  and  the  apparent  A>^-ef feet  less 
significant.  The  implication  of  this  is  that  by  increasing 
diffusion  time  in  a stirred-reservoir  experiment, 
the  error  due  to  the  -effect  is  minimised.  This 
follows  an  earlier  suggestion  made  by  Richards  (21).' 

If,  in  the  analysis  of  capillary  experiments 


Ji 


with  stirring  of  the  reservoir,  the  A 5-effect  is 
neglected,  the  actual  length  of  capillary  in  which 
diffusion-only  occurs  will  bo  shorter  than  the 
actual  length  of  the  capillary,  and  therefore,  the 
calculated  diffusion  coefficient  v;ill  appear  larger 
than  should  be  expected,' 

(b) . The  No-Stirrinr?  Noc’el 

In  order  for  there  not  to  be  a -effect", 
there  can  be  no  stirring  of  the  reservoir  solution 
outside  the  capillary.'  If  the  reservoir  material 
is  not  stirred,  hov/ever,'  the  capillary-reservoir 
interface  cannot  be  maintained  at  zero  difruscint 
concentration.'  Because  of  the  complexities  of 
solution  of  the  diffusion  equation  v;ith  a changing 
boundary  condition,'  experimental  results  based  on 
the  no-stirring  model  V7ere  not  reported  until  a 
computer-simulated  mathematical  model  was 
developed  (22),  Tlie  experiments  described  above, 
all  done  in  the  1950's  and  1960's,  v/ere  based  on 
the  stirring  model  of  capillary  diffusion,' 

B.  Existing  Djffusivity  Data 
1,'  Potassium  Chloride  (KCl) 


I 


(that  is,‘  infinite  dilution)  is  given  as  l,99xl0“5  cra^/sec 
at  (23). • The  minimum  value  is  given  as  1. '835x10-5 

at  0,2  molar  concentration;  the  diffusivity  then  increases 
with  increasing  concentration,'  At  1,0  molar  KCl,  the 
diffusion  coefficient  is  1,893x10-5  cm^/sec.'  C,Y,  Luk 
integi’atcd  the  dOita  described  above  to  obtain  an  average 
(integral)  diffusion  coefficient  for  1,0  molar  KCl 
diffusing  into  pure  u-ater  at  25°C  of  1,866x10“-^  cra^/sec,  (20) 


As  this  value  is  generally  accepted  as  valid,  it 
was  considered  to  be  a good  starting  point  to  do 
no-stirring  diffusion  experiments  v/ith  1.0  molar  KCl 


in  water  at  *.5^C  anu. 


cfciiu.  I/O  OOiViO^i'O  i/iiO  1‘OoUxuo  Wxi/ii  liiic 


value  listed  above. 


Zinc  Oxide  (ZnO) 

The  system  under  consideration  is  0.01  molar 
zinc  oxide  in  40/^  (v7t.')  potassium  hydroxide,  diffusing 
into  ^0%  KOK  at  25°C,’  This  electrolyte  mixture 
is  of  particular  interest  to  electrochemical  engineers, 
as  it  is  commonly  used  in  battery  systems,*  Data  have 
been  published  (24)  for  the  diffusion  coefficient  of 
O.'Ol  molar  ZnO  in  23.4?J  KOH  as  7.10x10”^  cm^/sec, 
along  with  ■values  for  other  concentrations  of  zincate 
ion,* 


f 


3.  Cupric  Sulfate  (CuSO/i) 

There  has  been  one  generally  accepted  "standard" 
set  of  values  for  diffusion  coefficients  for  cupric 
Eulf ate/sulfuric  acid  aqueous  solutions  (25),  These 
data,  have  been  around  for  about  forty  years,  and 
have  formed  the  basis  of  most  experiments  done 
in  that  time  for  v:hich  these  diffusivitics  are 
required.’  The  values  obtained  from  these  experiments 
by  Cole  and  Gordon,  shovm  in  Table  2,1,  were  obtained 
by  a diaphragm  cell  technique,  and  v;ere  run  at  18®C, 

A more  recent  study  by  Eversole  (26)  gives 

diffusing  in  vmter  at  25°C  (see  Table  2.2),  Both 

these  data  and  that  of  Cole  and  Gordon  (Table  2,1) 

indicate  a concentration-dependence  of  the  diffusion 

coefficient.  Yet,  as  discussed  later  in  section  III-A, 

nd 

the  solution  of  Pick's  2 Lavj  of  Diffusion  (the 
•diffusion  equation")  is  based  on  the  assumption  of 
constant  diffusivity.  It  may  be  argued,  though, 
that  the  change  in  diffusivity  vrith  respect  to 
concentration  change  is  small  enough  to  make  that  assumption 
reasonable  at  low  concentrations.  This  paradox  vrill 
be  discussed  at  the  end  of  this  report,  in  the  light 
of  the  data  in  Tables  2.1  and  2,2,  as  vfcll  as  the 
present  study. 


\ 

i 
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Table  2.1  - Data  of  Cole  and  Gordon  (25) 

for  Cu30ij,/H2S0i^,  solutions  diffusing 
into  H2S0ij.  at  18°C, 


Molarity 
of  CuSOjij, 

Molarity 
of  H2S0/(. 

X 106 

0.05 

0.00 

6.05 

0.10 

0.00 

5.85 

0.15 

0.00 

5.50 

0.05 

0.‘50 

5.95 

0.10 

0.50 

5.90 

0.15 

o.'50 

5.75 

0.05 

1.80 

5.00 

0.10 

1.80 

4.95 

0.15 

1 .80 

4.90 

1 


jfable  2._2  - Data  of  Eversole  ojb  al.  (26)  for 

CUSO/1/H2O  solutions  diffusing  into 
water  at  250C, 


Molarity  of  CuSOJ[. 

0,0028 

0.0056 

0.0084 

0.0112 

0.‘0308 

0,0980 

0.2100 


X 10^  cn^/^eo 

7.474 

7.159 
6.922 
6.726 

6.159 

5.644 

5.349 


— A 


f 


I 


[ 
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III,  Tlieorctjcal  BnclfrTOur^d 

A.  Solut'.ions  of  tho  Diffur^ion  Equation 
;!.*  Thft  Diffit.sion  ETuatlon 


We  are  concerned  with  dif fusion  of  electrolyte 
out  of  a capillary  into  a bulk  solution.  In  the  capillary 
is  the  electrolyte  under  study  in  a solution  of 
supporting  electrolyte  (vjator,'  acid,  or  base  solution) . 

The  bulk  liquid  (surrounding  the  capillary)  is  the 
supporting  electrolyte  only.  Temperature  is  assumed 
constant  in  and  around  the  capillary;  therefore, 
thermal  effects  arc  assumed  non-existent.  Diffusion 
i^ithin  the  capillary  is  assumed  to  be  one-dimensional. 

The  governing  diffusion  equation  is  knovm  as  Pick’s 
SecondLaw  Of  Diffusion: 


where, 


3c  _ ^ pZc 


(3.1) 


c = concentration  of  diffusant  in  capillary 
t = time 

X = distance  from  closed  end  of  capillary 
^ s diffusion  coefficient 


For  convenience  in  later  discussion,  the 
follov;ing  dimensionless  parameters  are  introduced: 


! 


! 


i 

i 


(|)  *=  c/cq  (cq  = initial  cone,  in  capillary) 

X s x/L  (L  = length  of  capillary) 

r=  t/i^ 

Equation  3.1  niay  now  be  rewritten  in  dimensionless  form: 


^6  _ 0 
OT 

The  initial  condition  is  defined  as : 


(3.2) 


(j)  = 1 at  , r=  0 (3.3) 


Boundary  conditions  v/ill  be  determined  by  v/hatever  is 
occurring  at  the  interface  of  the  bulk  solution  and  the 
capillary  mouth. 


It  has  been  demonstrated  by  Luk  (20)  that  a fluid 
velocity  of  2x10-^  cm. /sec  across  the  capillary  mouth  is 
sufficient  to  maintain  a zero-concentration  of  diffusant 
at  the  interface,  if  the  concentration  in  the  capillary 
is  low  (i.e.,  less  than  one  molar  for  CUSO4).  The  method 
of  calculating  the  bulk  fluid  velocity  (and  corresponding 
Reynold's  Number)  will  be  discussed  in  a later  chapter 
of  this  report.  It  may  be  assumed  that  this  minimum 
velocity  requirement  is  met,  and  that  there  is  no 
diffusing  electrolyte  outside  the  capillary  moutli. 


Thus,  the  boundary  conditions  for  equation  3.2  are: 


= 0 at  X=  0 , (3.*^) 

^ =0  at  >.=  1 , T>0  (3.5) 

Using  the  method  of  "separation  of  variables", 
equation  3.2,  v/ith  conditions  3,3,  3,^,  and  3.5 
may  be  solved  as : 


Note  that  as  dimensionless  time  (t)  increases,  the 
exponential  value  in  parenthesis  decreases.  For 
sufficiently  long  time,  all  but  the  first  tetm 
(n=0)  may  be  dropped.  At  L = 0.10,  for  example, 
neglecting  the  n=l  and  higher  terms  gives  a value 
of  P which  is  O.OIJ^  less  than  if  several  terms 
were  included,  AtT=  0,06,  however,  the  error  due 
to  neglecting  higher  order  terms  is  -3.^^,  and  is 
not  considered  negligible.  It  can  be  shown  that 
for  dimensionless  time  of  0,10  or  greater,  a 
single  term  equation  may  be  used  with  no 
appreciable  loss  of  accuracy.  It  is  possible, 
however,  to  solve  equation  3.2  by  other  methods, 
with  the  goal  of  finding  a solution  suitable  for 
use  with  small  dimensionless  time,  Cne  such  solution 


is  obtained  by  the  method  of  Laplace  Transforms,  as 
shown  in  equation  3,7; 


In  this  case,  it  can  be  shovm  that  the  entire  summation 
term  may  be  dropped  for  small  dimensionless  time  0.06). 
Thus,  for  capillary  diffusion  with  stirring  at  the 
capillary  mouth,  the  follov.’ing  tv;o  equations  are 
valid  in  the  specified  ranges  of  dimensionless  time: 

For  T ^0.10  , ()  = • EXP-Cn^rA)  (3.-8) 

Port  <^0.06  , (j)  = 1 - AT/rr)i-  (3.9) 

For  0,06<  T < 0,10,  several  terms  of  either  equation 
(3.6  or  3.7)  may  be  used.  These  equations  m^y  be 
rearranged  to  give  t as  a function  of  (J),  the  dimensionless 
concentration,'  Thus,  if  Cq  is  knovm,  capillary  length  (L) 
is  measured,  and  diffusion  time  (t)  is  specified,  it 
is  possible  to  compute  a diffusion  coefficient  from 
measurements  of  average  concentration  of  diffusant 
remaining  in  the  capillary  (c)  after  the  specified 
diffusion  time  period,' 

3.'  Calculation  of  A 2 


P' 


As  it  is  the  purpose  of  these  experiments  to 


determine  the  so-called  "a ^-effect"  in  stirred-rcservoir 

I 

experiments,  as  described  previously,  it  v/ill  be  useful 
to  have  a way  of  calculating  the  magnitude  of  If 

the  equations  of  the  previous  section  (3.6  through  3.9) 
are  used  to  calculate  a dimensionless  time  from  a 
measurement  of  concentration,  and  if  a diffusion 
coefficient  'is  knovm  (from  tabulated  values,  or 
from  other  experiments),  then  an  effective  diffusion 
length  (L')  may  be  calculated.'  L*  should  be  considered 
to  be  the  length  of  capillary  in  v;hich  diffusion  alone 
occurs.'  Since  the  real  length  of  the  capillary  is 
known  (L),'  and  the  effective  length  may  be  calculated 
(L*)  the  section  of  capillary  in  which  non-diffusive 
mass  transport  occurs  may  be  calculated  as: 

« L - L*  (3.10) 

Wie  results  of  no-stirring  experiments,  from  which 
a diffusivity  may  also  be  calculated^'  may  be  used 
to  calculate  in  the  stirred-reservoir  experiments. 

No-St trr in?  of  the  Bulk  Solution 

The  "a /-effect"  discussed  above  may  be  avoided 
by  not  stirring  the  bulk  fluid  around  the  capillary. 

If  there  is  no  stirring,  however,  the  zero-concentration 
boundary  condition  (equation  3.5)  is  not  valid.'  In  this 


caso,  it  is  necossary  to  consider  diffusion  outside 
, the  capillary  in  three  dimensions,  and  the  solution 

i of  equation  3.2  is  very  complicated,  A paper 

t 

entitled  "An  Improved  Model  for  Capillary  Diffusion 
with  No  External  Stirring"  (22)  discusses  computer 

I 

; modeling  of  this  situation,  and  resultant  numerical 

I solutions.  Figure  ^ is  a reproduction  of  the 

f plot  of  ^ vs,’  T from  the  paper  described  above,' 

I It  shov/s  both  the  analytic  solution  for  the 

j stirring  boundary  condition  (equation  3,6  or  3.7), 

: and  the  numerical  solution  for  no-stirring,'  (A 

« 

I third  curve  shov/s  an  analytic  solution  based  on 

I 

the  assumption  of  one-dimensional  diffusion  with 
, no  stirring,  but  this  appears  to  differ  greatly 

t 

I 

i from  the  numerical  solution  for  all  but  very  small 

I values  of  T ) .' 

I 

I . 

i 

j It  may  be  inferred  from  Figure  4 that,  in  the 

i limit  of  decreasing  dimensionless  time,  the  no-stir 

I numerical  solution  approaches  the  short-time 

( solution  for  stirring  (equation  3.9).’  By  assuming 

that  the  short-time  no-stirring  approximation  v;ill 

have  the  same  form  as  equation  3. ‘9,  and  plotting 


the  data  of  Figure  4 on  a log-log  scale,  the  following 


Average  capillary  concentration  (Oave^ 
vs.  dimensionless  time  (t) . (22) 

• • No  stir,  three  din,’,  computed 

— — — - Stirring,  one  dim, , analytic 
— No  stir,  one  din,'  , analytic 


empirical  expression  is  obtained; 


s 1 - ; 0<r<0.3  (3.11) 

This  equation  coincides  v/ith  the  no-stirring  numerical 
solution  to  less  than  1%  deviation  for  the  range  of 
dimensionless  time  listed  above,  T}ius,  short-time 
experiments  may  be  run  with  no  stirring  of  the  bulk 
solution,  and  diffusion  coefficients  may  be  calculated,* 
(It  should  be  noted  that  a similar  procedure  to  the 
above  could  be  used  to  obtain  a long-time  approximation 
for  no-stirring  experiments,-  also.*  All  experiments 
described  in  this  i-*eport,  however,  vjere  run  in 
the  bxjwL— I/iu«#  (lefliieu.  above,’) 


IV.  ExHorlrontal  Apparatus  and  Procedures 
• A.  CarlllarlGS 

T»Jo  sets  of  luclte  capillaries  v;ere  used;  one 
set  for  the  no-stirring  experiments,  and  the  other  for 
those  experiments  during  which  the  reservoir  was 

« 

I stirred.  For  the  no-stirring  mode,  the  capillaries  used 

I were  the  same  as  described  by  C,  Y.  Luk  in  his  research 

t thesis  (20),  Each  of  the  four  capillaries  looked  like 

t the  one  shovm  in  Figure  5i  except  for  differences  in 

internal  length  (shaded  area  in  figure,  notation  "L"), 
IVo  of  the  no-stir  capillaries  had  internal  length 
of  one  inch  (2,5^  cm.)  and  the  other  two  were  of 
1.3^5  Inch  (3.49  cm.)  internal  length.  The  internal 
diameter  of  each  capillary  was  l/l6  inch  (0,158  cm,). 

The  flat  plate  and  collar  at  the  mouth  of  the  capillary 
(shown  in  the  figure)  act  to. minimize  convection  caused 
by  diffusant  which  would  fall  down  the  side  of  the 
i capillary  if  no  collar  was  present. 

The  need  for  the  collar  is  demonstrated  in 
Figures  6a  and  6b.^*''\n  Figure  6a,  a capillary  is 
. initially  filled  with  concentrated  KMnO^  and  is 
placed  in  a water-filled  reservoir.  Even  without 
external  stirring,  the  denser-than-water  KMnO^ 


L 


Bolution  falls  over  the  side  of  the  capillary, 
creating  an  undesirable  stirring  effect.  In 
Figure  6b,  the  same  capillary  is  shown,  this  time 
with  a collar-and -barrier  arrangement,  held  in 
place  v/ith  black  wax.  Note  that  nov/  the  diffusing 
KMnOii  solution  remains  in  the  vicinity  of  the 

I 

capillary  mouth,  as  described  by  the  no-stirring 
boundary  conditions,  (In  the  figure,  a small 
amount  of  KMn04  seems  to  be  falling  avray  from  the 
capillary.  This  is  due  to  a poor  seal  vjith  the 
wax.  New  capillaries  were  made  v;ith  the  collar 
an  integral  part  of  the  whole  capillary,  and 
these  new  capillaries,  as  shovm  in  Figure  5»  v;ere 
the  ones  used  in  the  no-stirring  experiments.) 

For  the  stirred-reservoir  experiments,  four 
new  capillaries  were  made,  each  1,0  cm.  in  internal 
length  and  0,1  cm.  internal  diameter.  These 
capillaries  look  like  the  one  shown  in  Figure  6a, 
There  is  no  collar  on  the  capillary  in  this  series 
of  experiments,  because  forced  convection  is  the 
desired  result  of  stirring  the  reservoir,  unlike 
the  no-stirring  model,  in  which  all  convection  should 
be  eliminated. 


Drawing  (rot  to  scale)  of  capillary 
used. in  no-stirring  experiments. 
Diffusant  is  in  shaded  area. 


Flppjre  6a  - KHnO^-filled  capillary  ' Figure  6b  - Same  capillary  as 
in  water  reservoir,’  Note  that  at  left  except  that  anti- 

KMnOii,  falls  due  to  convection,*  convection  ‘carrier  is  in  place 
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Figure  6a  - KMnO^-filled  capillary  ' Fip:ur^.6;b  - Same  capillary  as 
in  water  reservoir.'  Note  that  at  left  except  that  anti- 

KMnOi^  falls  due  to  convection.'  convection  carrier  is  in  place 


B.  Floi7-Cbanncl  / Rogcrvotr 


A rectangular-shapsd  lucite  box  was  used  to 
hold  the  capillaries  and  surrounding  reservoir 
material,'  The  box  v:as  designed  so  that  it  could  be 
used  as  a flow-channel  for  stirred-reservoir 
experiments  or  as  a static  reservoir  for  the 
no-stirring  experiments,'  Holes  were  made  in  the 
base  of  the  reservoir  box  so  that  up  to  six 
capillaries  at  one  time  could  sit  in  the  reservoir. 

The  outside  dimensions  of  the  capillaries  and 
the  locations  of  the  six  holes  were  made  so  as 
to  to  insure  that  the  top  of  each  capinary  v;onlrt 
be  at  the  same  height  and  distance  from  sidev:all3 
as  all  of  the  others.  This  arrangement  ass*. res 
that  in  the  stirred-reservoir  experiments,  the 
fluid  velocity  passing  the  capillaries  vrould  be  equal 
at  each  capillary-reservoir  Interface.  As  shovm  in 
the  figure,  the  flov;  pattern  is  defined  by  the  geometry 
of  the  box,  v;ith  incoming  reservoir  fluid  entering 
at  the  opening  labelled  "A"  and  leaving  at  "B"  (see 
Figure  7).'  A thermal  regulator  maintains  constant 
temperature,  and  a constant-head  tank  maintains 
constant  flow  rate.'  The  overall  flow  chart  is 
sho\m  In  Figure  7 on  the  next  page. 


The  no-stirring  experiments  do  not  require  such 
an  elaborate  flov/  system  as  described  above.  All  that 
is  required  is  a box  which  holds  the  capillaries  in  a 
static  reservoir,  all  at  uniform  temperature. 

In  order  to  insure  a uniform,  non-convective 
start-up  of  the  no-stirring  experiments,  nylon  screv/s 
were  installed  in  the  top  of  the  reservoir  box,  located 
above  each  capillary.  When  an  experiment  v;as  to  begin, 
the  diffusant-filled  capillary  v;as  placed  in  the 
partially-filled  reservoir  box.'  (The  level  of  the 
reservoir  had  to  be  belovf  the  capillary  mouth,  so  that 
no  splashing  of  the  reservoir  fluid  would  remove 
capillary  material.)  The  top  v;as  firmly  placed  atop 
the  reservoir  box,  and  the  nylon  scrov;s  were  tightened, 
forming  a seal  at  the  capillary-reservoir  interface. 

At  this  time,  the  reservoir  could  be  filled  to  the 
top,  and  the  experiment  v;ould  begin  vrhen  the  nylon 
screws  v^ere  removed,  thus  exposing  the  capillaries  to 
the  reservoir, 

C.  Temoerature  Control 
1.*  Stirrin.fr  ftcoerlments 

As  shovm  in  Figure  7 (p.  32 ) , the  temperature  of 
the  reservoir  vrater  is  maintained  at  a predetermined 


constant  value  by  a feedback  control  system.  In  a 
thirty-liter  holding  tank,  a temperature  regulator 
maintains  constant  temperature,  with  a stirring 
blade  in  the  tank  keeping  the  vrater  moving.  A 
circulating  pump  sends  the  heated  v;ater  to  a two-liter 
beaker  which  acts  as  a constant-head  tank,  sending 
water  at  constant  pressure  (i.e.,  constant  flow 
rate)  into  the  reservoir  box.  Out  of  the  reservoir 
box,  the  water  is  returned  to  the  heating  tank. 

Temperature  is  monitored  with  a copper  and 
constantan  (in  thin  glass  tubing)  thermocouple, 

InsertRd  into  the  reservoir,  oonnected  to  a 
laillivoltmetcr  and  chart  recorder,'  Two  thermo- 
couples t-jere  connected  in  series,  v/ith  another 
pair  of  thermocouples  at  an  ice-and-v;ater  (OOC) 
reference,  all  read  on  a Hewlett-Packard  digital 
voltmeter,  allowing  readings  accurate  to  0,5°C, 

2,*  No-St irrinrr  Pxneriments 

During  those  diffusion  experiments  in  vjhich  fluid 
does  not  flovr  through  the  reservoir  chamber  (that  is, 
no  stii‘ring),  the  temperature  of  the  static  reservoir 
is  kept  constant  by  holding  the  entire  reservoir 
box  in  a pool  of  constant  temperature  i-iater.  Thus 


tho  circulating  system  vihich  maintains  constant 
temperature  in  the  reservoir  bath  for  stirrcd-recervoir 
experiments  is  also  used  for  tho  no-stirring  series, 
except  that  tho  heated  v.viter  surrounds  the  box, 
rather  than  flowing  through  it,‘  Again,  temperature 
is  monitored  throughout  the  diffusion  experiments.’ 

D.  Analytical  Toehn^aue 

When  a diffusion  experiment  is  finished,  each 
capillary  must  be  removed  from  the  reservoir  and  its 
contents  analysed  to  determine  the  average  concentration 
••  of  diffusant  remaining  after  the  allotted  diffusion  time. 
The  analytical  device  used  in  this  determiriation  was 
a Techtron  model  AA-120  Atomic  Absorption  Spectro- 
photometer (A.A.S.),  which  utilizes  the  principle  of 
optical  absorption  of  a specific  x/avelength  corresponding 
to  a specific  element,  giving  a measurement  of  ionic 
concentration  read  on  the  A.A.S.  in  the  form  of  "absorbance 
imits"  (on  an  arbitrary  sero-to-hundred  scale).  The 
rocoramended  range  of  concentrations  to  be  used  with 
the  A.A.S.  , for  the  elements  of  interest  in  these 
experiments,  \7as  on  the  order  of  10”^  molar  ion 
concentration.  In  order  to  calibrate  the  A.A.S.,  a 
scries  of  standard  solutions,  across  a concentration 


i 


range  of  10  ^ to  10  molar,  \ms  made  up  for  each  of 
the  elements  to  be  analyzed  (Zn,  Cu,  K) , 

In  order  to  analyze  the  capillary  contents  after 
diffusion,  it  is  necessary  to  reduce  the  ion  concentration 
to  the  working  range  of  the  A.A.S.  This  is  accomplished 
by  d^uting  the  capillary  contents  into  a Jaiovm 
volume  of  distilled  vmxter.  V/hen  converting  from  the 
■absorbance  units"  readout  on  the  A.A.S.  to  a molar 
concentration,  this  dilutirn  factor  is  accounted  for. 

E.  Calculation  of  Canillarv  Volumes 

Capillary  volumes  v:ere  determined  by  running  a 

■zero-time"  experiment;  that  is,  by  filling  each 

capillary  v^ith  diffucant  of  known  ion  concentration 

and  immediately  emptying  and  analyzing  the  contents. 

If  no  diffusion  vras  allovfed  to  occur,  then  the 

apparent  concentration  measui’ed  on  the  A.A.S.  should 

agree  with  the  laiowm  concentration  of  diffusant 

0 

initially  in  the  capillary.  Thus  the  only  unknovm 
value  is  the  capillary  volume,  v;hich  may  then  be  i 

f 

calculated  directly.  i 

P.  Calculation  of  r?c\Tiold*s  Number  (Stirring  Experiments) 

1 
1 

i 


i 

1 


Fluid  velocities  v;ere  calculated  by  assuming 


laminar  flow  in  the  flow  channel,  with  kno;m  ceoractry,' 
The  velocity  profile  in  the  horizontal  direction  is 
given  by  equation  4-.1  belov;,  with  Vq,  the  centerline 
velocity,  described  by  equation  .• 

V(y,z)  = Vc[^l  - (2y/w)5  (^.1) 

Vo  = - (z /D)5’  (4.2) 

where,’  V©  = centerline  velocity,  cm/sec 

Vt  * surface  velocity  above  centerline 
w « channel  widih,  era 
D e water  depth  in  channel,  cm 
y « horizontal  distance  from  centerline,  era 
z s!  vertical  distance  '^rom  cen'^orl^r'^ 

The  average  velocity  in  the  flov;  channel,  V,  is 
based  on  and  the  assumption  of  laminar  flow,  as 
given  by  equation  4.3.  This  average  velocity  may  also 
be  expressed  as  the  ratio  of  volumetric  flov/  through 
the  channel  (Q)  to  the  cross-sectional  area  of  flow  (A) 

Yt  ■=  9/4V  = 9/4*Q/A  (4.3) 

The  follov/ing  values  have  been  measured  in  the 
flow  channel:  The  capillaries  are  located  at 
y *®  +0.7  »’2  = 0.5.  The  dimensions  of  the  channel 
are  such  that:  v/  * 4.6  , D = 1,8  .•  Insertion  of 


these  values  into  equations  k,l,  4.2,  and  4.3  gives 
rise  to  the  following  expression  for  V(y,z)  as  a 
function  of  flov/  rate  (Q)  for  this  Imovm  geometry: 

V(y,z)  = 0.2275  Q (4.4) 


Thus,  the  Reynold's  Number  at  the  capillary  mouth 
may  be  calculated  from  a measured  flov/  rate,  according 
to  the  following  formula: 


Re 


- 0.2^75  d /O  0. 


(4.5) 


in  which: 


d = capillary  inside  diameter,  cm 
= density  of  reservoir  fluid,  gm/cm3 


* V » W ate*  A ^ V it  A *<««  V Lr  ^ ^ •&«««'  / O ^ 


)X  - reservoir  fluid  viscosity,  poise 


G.  Experimental  Frocorlure 


1.*  No-St5.rrin'r  S::rer iment 

a. '  Pill  reservoir  partly  (about  75/'i  full) 

b. ’  Place  reservoir  box  in  constant- temperature 

bath;  wait  until  internal  temperature  reaches 
steady-state.' 

c. '  Pill  the  capillaries  with  diffusant,  using 

a sterile  hypodermic  needle  for  each. 

d.  Insert  capillaries  into  reservoir  box 

and  close  box,'  Seal  tops  with  nylon  screv/s. 


e, '  Fill  reservoir  completely;  v/ait  again  for 
. temperature  to  stabilize. 

f, '  To  initiate  diffusion,  slov/ly  turn  the 

nylon  screv/s  until  each  is  as  high  as  it  can 
go  without  coming  out  of  the  box,  Tlie 
capillary-reservoir  interface  is  nov/  free 
of  obstruction  and  diffusion  v;ill  proceed 
imtil  the  screv;s  are  tightened  down, 

g, '  After  diffusion  is  stopped  (nylon  screws 

back  in  place)  and  time  is  noted,  empty 
the  reservoir  box  of  fluid,  by  removing 
the  box  from  the  constant-temperature 
bath,  and  opening  one  of  the  holes  in  the 
side  of  the  box,  allowing  reservoir  fluid 
to  drain  out. 

h, ’  V/hen  the  level  of  the  reservoir  material  is 

well  belov;  that  of  the  capillary  tops, 
remove  the  top  of  the  reservoir  box  and 
then  remove  each  of  the  capillaries  from 
the  reservoir. 

1.'  Using  a fresh,  sterile  hypodermic  needle 
(use  a fresh  one  for  each  capillary) , remove 
the  capillary's  contents  and  place  in  a 
volumetric  flask.  Fill  the  volumetric 


i 


flaslc  up  to  the  nark,  first  with  the  washings 

. of  the  capillary  (wash  several  times  with 

water)  and  then  with  v/atcr  for  the  remainder. 

The  capillar3'^  contents  have  nov;  been  removed 

and  diluted  to  a Imovm  volume,  as  described 

earlier  in  this  section  (see  p.  36), 

« 

J,‘  Proceed  v/ith  analysis  on  the  A.A.S.,  to 
obtain  a measured  value  of  average  concen- 
tration of  ion  remaining  in  the  capillary 
after  diffusion.  This  value,  divided  by 
the  initial  ion  concentration,  is  and 
may  be  used  in  the  appropriate  equation 
found  on  pp,  21-22  to  find  a diffusion 
coefficient. 

2,‘  Stirred-Reservo'^.r  Ercneriment 

a, '  Fill  reservoir  and  start  flow  system. 

Allow  temperature  to  stabilize. 

b, '  Periodically  check  flow  rate  both  before 

and  during  experiment, 
o.  Pill  capillaries  with  diffusant,  using  a 
sterile  hypodermic  needle. 

d.  Insert  capillaries  in  reservoir,  noting  time. 

e. ‘  After  desired  diffusion  time  is  passed,  remove 

capillaries  from  flow  channel. 

f. *  Continue  steps  "i”  and  "j"  from  above. 


V.  Exnerlmontal  Results 

A.  No-St irrin.r  E::r?qrlmontn 

1,*  Potassium  Chlorj.de  In  V/ater 

The  initial  no-stirring  diffusion  experiment 
involved  l.OM  KCl  in  H2O  diffusing  into  H2O  in  a 
static  reservoir  at  a temperature  of  25°C,  Using 
four  capillaries  v/ith  the  anti-convection  collars, 
the  results  listed  in  Table  5.1  were  obtained.  As 
shown,  the  average  diffusion  coefficient  of  the  four 
data  points  is:  1.843  x 10“5  cm^/sec.  The 
standard  deviation  is  0.183.x  10“5 

ZJ  Zinc  Oxide  in  Potassium  Hydroxide 

A series  of  diffusion  experiments  ;vere  run  in  which 
0,01  molar  ZnO  in  405^  (by  vieight)  KOH  diffused  into 
UO^  KOH  at  25®C,  in  a non-stirred  reservoir.  The 
results  of  these  experiments,  two  series  of  four 
capillaries,  are  listed  in  Table  5.2,  showing  an 
average  diffusivity  of  6.45  x 10"^  cm^/sec,  and 
a standard  deviati  .>n  of  2,42  x 10”^,  The  tabulated 
data  is  listed  in  tv/o  groups,  according  to  the 
length  of  the  capillary. 


Cupric  Sulfate  in  V/atcr 


Tabic  5.3  shows  the  results  of  experiments 
in  which  0.05  molar  CuSOj^  in  aqueous  solution 
diffusing  into  static  HgO  at  25°C,  Again,  the 
data  are  divided  into  tv;o  groups  according  to 
capillary  length.'  As  indicated,  the  average 
diffusion  coefficient  for  the  eight  tabulated 
values  is  5.560  x 10~^  cm^/sec,  v;ith  a standard 
deviation  of  2,066  x 10“^. 

4.'  Cuoric  Sulfate  in  Sulfuric  Acid 

The  CuS0|^/H2S0j^  system  was  studied  over  a 
range  of  cupric  ion  concentrations  and  at  tv;o 
temperatiu'es.  Cupric  sulfate,  in  a solution 
of  half-molar  sulfuric  acid,  diffused  for 
24  hours  in  a static  reservoir  of  0,5M  K2S0^, 

At  25®C,  cupric  concentrations  of  0.005?!,  0.02M, 
0,‘05M,  and  0.1  OM  were  studied.  At  35.5*^C,  the 
CuSOji^  concentrations  v.’cre  0.005M,  O.05H,  and  O.IOM. 
For  each  temperature  and  concentration.  Table  5.4 
lists  the  number  of  experimental  data  points  ("N"), 
the  concentration,  temperature,  average  diffusion 
coefficient,  and  standard  deviation. 


1^.  Stirrcd-Hosorvotr  S:vT?oringnfcs 

In  the  experiments  in  which  the  reservoir  was 
stirred,  an  aqueous  solution  of  0.05  molar  CuSOji^ 
at  26°C  initially  filled  the  capillaries,  v;hich 
sat  in  a reservoir  of  flowing  H^O,  Two  sets  of 
experiments  iiere  run:  one  at  a flov;  rate  of 
0,817  ral/sec  “ 2,03)  and  the  other  at 

2,'915  ral/sec  (Ng^  = 7.55).  In  each  series, 
several  diffusion  runs  were  made,  v;ith  the 
diffusion  times  varying  from  tv:o  to  twenty-four 
hours 

Table  5.  4 e«  n 

showing  for  each  flov^  rate  and  diffusion  time,  the 
average  diffusion  coefficient,  the  standard  deviation, 
and  the  number  of  data  points. 


TABLE  5.1  - Experimental  Results  for  1.0  Molar  KCl 
In  HgO  diffusing  into  H2O. 


L.  cm. 

t.  Mrs.  ^aye 

r 

105 

2.54 

'-’O 

24.0  0.481 

0.2546 

1.901 

2.54 

24.0  0.484 

0.2515 

1.878 

3.49 

24.0  0.644 

0.1122 

1.582 

3.49 

24.0  0.605 

0.1427 

2.011 

. 

Reservoir  not  stirred. 

Capillary  ^ntarnal  d^a^1(^^.aT':  0. 

1 58  em . 

Reservoir  temperature: 

25.0  + 0 

.5  ®c. 

^avg.  ” X 10“^  cm^/sec 

Standard  Deviation  = 0.183  x 10-'5  cir.2/sec 

95%  confidence  limits:  1.73  < J&'  <1.95  (x  10-5) 


TABLE  5.2  - Experimental  Results  for  0,01  Molar  ZnO 
In  40^  (v/t.)  KOH  diffusing  into  KOH, 


X#«  OTTl  • 

t.  hrs. 

Cave 

Co 

r 

106 

2.54 

23.0 

0.6764 

0.0986 

7.68 

2/54 

23.0 

0.7089 

0.0811 

6.32 

2.54 

24.5 

0.8877 

0.0141 

1.04 

2.54 

24.5 

0.7089 

0.0811 

5.93 

3.49 

23.0 

0.7361 

0.0677 

9.95 

3.49 

23.0 

0.7475 

0.0624 

9.17 

3.49 

24.5 

0.7655 

0.0545 

7.68 

3.49 

24.5 

0.8373 

0.0278 

3.84 

Reservoir 

not  stirred. 

Capillary 

internal  diameter : 0 

.158  cm. 

• 

Reservoir 

temperature: 

25.0  ± 

0.5  oc. 

^avg.  “ ^ cm2 /sec 

Standard  deviation  = 2,90  x 10"^  cm^/sec 

95^  confidence  limits:  ^.03  ^^8,8?  (x  10“^) 


TABLE  5«3  - Experimental  Results  for  0,05  Molar  CuSOij, 
in  H2O  Diffusing  into  H2O. 


L.  cm. 

t.  hrs. 

Cave 

r 

^x  10^ 

2.54 

22.5 

Cq 

0.7695 

0.0420 

3.344 

2.54 

22.5 

0.6305 

0.1258 

10.020 

2.54 

24.0 

0.7025 

0.0844 

6.303 

2.54 

24.0 

0.7021 

0.0846 

6.321 

3.49 

22.5 

0.8328 

0.0292 

4.396 

3.49 

22.5 

0.8327 

0.0293 

4.398 

3.49 

24.0 

0.8105 

0.0368 

5.191 

3.^9 

24.0 

0.B244 

0.0319 

4.510 

Reservoir 

not  stirred. 

Capillary 

internal  diameter: 

0,158  cm. 

Reservoir  temperature: 

25.0  ; 

t 0.5  oc. 

^avg.-  " 5. 5^50  X 16“^  cmVsec 
Standard  deviation  = 2,066  x 10”^ 

cm^/sec 

95%  confidence  limits:  4.08  ^^^7.04  (x  10"^) 


TABLE  5.^  - Experimental  Results  for  CuSO^  in  0.5*1  ^£^^4 
Dif fusing  into  0,5  Molar  H2SOj^  , 

JL  (molar)  T.  a vr-  x 10^  cm2 /sec  std.  dev. 


10 

0.005 

25.0 

7.24  0.29 

1 

7 

0.02 

25.0 

5.12  0.29 

1 

9 

0.05 

25.0 

6.04  0.13 

12 

0,10 

25.0 

4.32  0.10 

1 

6 

0.005 

35.5 

7.63  1.13 

■ 

9 

0.05 

35.5 

6.92  ' 0.36 

8 

0.02 

35.5 

6.92  0.13 

a 

Reservoir 

not  stirred. 

Capillary 

internal  diameter: 

0.158  cm. 

Capillary 

lengths:  2.54  cm. 

1 3.49  cm.  (mixed) 

Reservoir 

temperature:  as  above,  + 0,5°C, 

, 

N s=  number  of  data  points  represented 

Cq  “ initial  concentration  of  CuS0|^  in  capillary 


TABLE  5»5  - Experimental  Recults  for  0.05  Molar  CuSOj^ 
in  V/ater  Diffusing  into  ?lov;inr  VJater. 

Ji.  ^Re  Time.  Mrs,  x 10^  cn^/sec  std.  dev 


8 

2.0 

2 

7.877 

0.47 

7 

2.0 

3 

6.317 

1.07 

2.0* 

10 

6.453 

0.13 

8 

2,0 

18 

5.763 

0.62 

7 

2.0 

24 

6.088 

0.07 

4 

7.6 

2 

6.699 

0.67 

4 

7.6 

3 

7.866 

0.31 

4 

7.6 

9 

4.746 

0.47 

4 

7.6 

18 

5.488 

0.05 

4 

7.6 

24 

5.520 

0.03 

Reservoir 

stirred 

, Reynold's  Number  given 

above , 

Capillary  intci'nal 

diameter:  0,100  cm. 

Capillary  length: 

1,0  cm. 

Reservoir 

temperature:  26.0  ± 1,0  °C. 

N « number 

■ of  data 

points  represented 

VI.  Disciifssjon  of  Data 


A,  No-Stirr5r>r:  R>;T?3r1.mcTits 

1.'  Potassium  Chloridb  Ir)  VJater 

The  data  in  Table  5.1  represent  one  diffusion 
experiment  in  v;hich  four  capillaries,  each  filled 
with  a one-raolar  solution  of  KCl  in  distilled 
water,  sat  in  a 25°  water  reservoir  for  24-  hours, 
with  no  external  stirring  of  the  reservoir.  The  mean 
diffusion  coefficient  of  the  four  capillaries, 
found  by  summing  the  four  individual  values  and 
dividing  by  the  number  of  data  points  (4-)  is: 

0 = 1 ,84-3  X 10"5  cm^/sec,'  The  standard  deviation 
of  the  data  is  0.183.  According  to  a double-tailed 
t-test  of  significance  at  the  0.95  confidence  level, 
for  four  data  points  (see  Appendix  C for  examples 
of  this  type  of  calculation) , any  value  of  0 such 
that  1,73^^$  1.95  is  considered  in  good  agreement 
with  the  experimental  data  of  Table  5.1. 

As  discussed  in  section  II-B,  the  accepted  value 
for  the  integral  diffusivity  of  IM  KCl  in  water  at 
25°C,  is  1.866  X 10"^  cra^/sec,  which  is  in  excellent 
agreement  with  the  above.  In  addition,  C.Y,  Luk  (20) 
found  values  of  1.86^.^^  1.89  in  his  experiments  with 


the  same  system.  The  95^  confidence  interval  of 
Luk*s  data  appears  narrov;er  than  that  given  in  the 
paragraph  above  for  this  study,  because  Luk  had 
36  experimental  data  points  to  evaluate,  whereas 
only  four  data  points  are  included  in  this  study. 

The  .potassium  chloride  experiment  was  intended 

as  an  evaluation  of  this  technique,  to  insure  that 

all  equipment  v/as  in  v/orking  order,  and  that  the  revised 

analytical  technique  (sec  pp,  23-26)  was  Valid, 

As  the  results  of  this  experiment  v/ere  in  such 

good  agreement  v^ith  previous  experimental  work 

as  well  as  the  conerallv  "acc anted " value,  it  i-zas 
^ • • 

decided  to  continue  v/ith  this  method  for  zincate 
and  cupric  ion  diffusivities. 

2.  Zinc  Oxide  in  Fotassium  Hydroxide 

The  data  in  Table  5.2  represent  tv;o  diffusion 
experiments  capillaries  in  each  run)  in  v;hich  a 
solution  of.  O.OIH  ZnO  in  KOH  diffused  into  a 
static  reservoir  of  ^■0^  KOH,  The  mean  diffusion 
coefficient  and  the  95)^  confidence  limits  of  the 
data  are : 

^avg,  “ X 10“^  cmVsec 

Note  tliat  one  point  in  Table  5.2  lies  outside  of 


{ 


this  limit  (at  ^ - 1,04  x 10“^),  Further,  as  this 
point  is  even  outside  the  confidence  interval 
of  the  data,  it  may  be  ncglcctcdd,  and  a nev/  mean 
value  calculated.  The  new  mean  (of  the  seven 
remaining  points)  and  corresponding  95^  confidence 
limits  are: 

= 7.22  + 1.90  X 10-^  cmVsec 

As  the  range  of  the  scatter  is  narrov/er  when  the 
extraneous  point  is  not  included,  the  standard 
deviation  is  lov/er  and  the  limits  of  confidence 
are  reduced. 

This  diffusion  coefficient  of  O.OIM  ZnO  in 
40/^  KOH  may  be  compared  with  the  data  of  Lu  (24) 

(as  described  in  section  II-3-2),  Ke  studied  the 
system  of  O.OIM  ZnO  in  23, 45^  KOH  end  found  a 
diffusion  coefficient  of  7.10  x 10“^  cm^/sec  (no 
range  of  error  reported) . According  to  a double- 
tailed t-test  of  significance,  within  95^  confidence 
limits,  there  is  no  difference  betv;een  Lu*s  value 
and  the  average  diffusion  coefficient  found  in  this 
study.  The  Stokes-Einstein  Equation  (bolov;) , however. 
Implies  that  the  two  diffusivities  should  be  different: 
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yj'-  RT  1 

^ “ K • 67fr  u 

where,  V - diffusion  coeff iceint,  om2/sec 
R = gas  constant,  cal/gm.mole-OK 
T « absolute  temperature, 
r = ionic  radius,  S 
u = viscosity,  of  medium,  cp, 

K = 6,02  X 10^3  (Avogadro’s  Number) 


(6.1) 


The  textbook  "Alkaline  Storage  Batteries"  (28)  gives 
the  viscosity  of  40^  KOK  solutions  as  3.5  ccntipoise 
and  of  23.4/'J  solutions  as  1.8  centipoise.  If  this 
is  so,  and  all  other  factors  in  equation  6,1  are 
constant,  then  the  expected  diffusivity  of  zincate 
in  40^  KOH  should  bo  roughly  half  of  that  of  the 
same  zincate  concentration  in  23.4?J  KOH,  As  it  turns 
out,  there  is  no  significant  difference  between  the 
two  experimentally-determined  values  within  95^ 
confidence  limits.  Possible  explanations  of  this 
are;  ionic  effects  (caused  by  increase  in  pH  corresponding 
to  higher  concentration  of  KOH,  possibly  affecting  the 
value  of  "r"  in  equation  6,1);  and  temperature 
differences  (Lu  states  in  his  report  only  that  the 
experiments  were  run  at  "room  temperature",  which 
could  mean  anything  in  the  range  of  20-2 5OC.), 


yj  - 


C.Y.  Luk  (20)  also  studied  the  O.OIM  ZnO  in 
KOH  system  and  found  a diffusion  coefficient, 
expressend  in  terms  of  his  confidence  interval 
as  6.00  X 6,68  x 10""^,  which  (according 

to  the  double-tailed  t-test  at  P = 0.05)  is  not 
different  from  the  data  obtained  in  this  study. 

The  data  for  the  ZnO/KOH  system  v.'ere  also 
studied  in  order  to  evaluate  any  relationship 
between  capillary  length  rnd  apparent  diffusion 
coefficient , in  the  no-stirring  experiment.  A 
double-tailed  t-test  of  significance  at  P = 0.05 

showed  no  di  ■rfe'Pf>nee  betw#»en  fVip  f.wn  rrt»r>imes  rtf' 

data  (at  L = 3.^9  and  L = 2.54  cm.),  indicating 
no  variation  in  ^ of  this  series  of  experiments 
due  to  capillary  length. 


Cuoric  Sulfate  in  Water 


The  data  in  Table  5.3  were  evaluated  in  the 
same  manner  as  described  above  for  KCl  and  ZnO/KOH. 
The  mean  diffusion  coefficient  v;as  found  to  be: 

5.56  ± 1,43  X 10"^  cm^/sec 

®vg« 

(expressed  with  95;^  confidence  limits).  One  point 
In  the  table  (at  ^ = 10.02  x 10"^)  is  outside  the 
9S>^  confidence  interval,  and  may  be  discarded.  This 


» 


I 

r 


lowers  the  mean  value  of  % and  reduces  the  scatter: 


^Avtr  ~ ^♦'•923  i 0.968  X 10“^  cm^/sec 

As  explained  in  section  II-S-3,  the  data  of 
Eversolc  ^ al.  (26)  represent  point  volues  of 
diffusivities  across  a large  range  of  concentrations 
(0  to  0.35  molar).  These  values,  often  called  the 
"differential  diffusion  coefficients"  must  be 
integrated  as  a function  of  concentration  to  ob'tain 
the  integral  diffusivity,  vjhich  is  the  type  of 
diffusion  coefficient  obtained  by  the  capillary 
method,  Eversole  et  al.  suggest  a relationship  of 
tue  form:  ^ — a r bC(y  , Tov  less  Lljaii  ,35**  CuoOi;.  1x1 


water.  In  the  range  of  concentrations  of  interest  in 
this  study  (Cq  O.O^H),  the  diffusivity  may  be 
approximated  as  a linear  function  of  concentration 
with  no  worse  than  2^  deviation  from  Eversole 's 
data.  The  data  of  Eversole  graphically 

Integrated  according  to  the  follov;ing  formula: 

^int.  “ ^ (6.2) 

where  ~ integral  diffusivity  at  C = Cji 

C5,  « concentration  of  CuS04  (C^  = 0.05  H) 

- point-value  of  from  tabulated  values 
by  Eversole  et  al_, , range  to  Ci+j , 


This  graphical  integration  of  the  data  from  = 0 
to  Ci  - ~ 0*05  molar  gives  a value  for  the 

integral  diffusivity  of  cupric  sulfate  in  v.’ater 
at  250  of  5«85  X 10“^  cm^/sec.  According  to  a 
double-tailed  t-tcst  of  significance  at  P = O.05, 
this  value  is  in  good  agreement  v;ith  the  experi- 
mentally-obtained  data  shovm  in  Table  5.3» 

4.  Cunric  St^lfate  in  Sulfuric  Acid 

Each  value  of  ^ lifted  in  Table  5«4  is  an 
average  of  a number  of  individual  data  points. 

There  are  two  groups  of  numbers;  four  values 

c*  uv«*i/c*jL  ^ ^ 

and  three  average  values  representing  a total  of 
23  data  points  taken  at  35.5°C.  Each  mean  value 
in  the  table  is  listed  with  the  standard  deviation. 
Bach  group  of  values  at  one  temperature  viere 
evaluated  to  determine  the  effect  of  concentration 
on  the  diffusion  coefficient.  A linear  regression 
of  the  38  points  at  25°C,  with  X = concentration 
and  Y = diffusivity  x 10^,  produced  the  following 
linear  equation:  Y = -12.812X  + 6.017.  However, 
an  P-test  of  significance  at  the  0.95  confidence 
level  with  ^6  degrees  of  freedom  indicates  no 
correlation  of  the  data.  Within  the  95/<>  limits  as 


described,  slope  of  the  linear  regression  could  be 
anything  from  -28,504  to  +2,830,  As  that  range 
includes  the  slope  of  zero,  v;hich  corresponds  to 
no  concentration  effect,  no  statement  may  be  made 
to  the  effect  that  diffusivity  is  a linear  function 
of  concentration,  based  on  these  experiments. 

Similarly,  for  the  23  points  at  35»5°C, 
a linear  regression  gives  a line  of  the  form: 

Y « -4.157X  + 7.205.  Again,  within  the  limits 
of  the  95/^  confidence  level,  the  upper  and  lo:;er 
bounds  of  the  slopes  of  the  regression  line  are, 
respectively,  16,732  and  -25.046,  crossing  the 
X-axis. 

Following  the  suggestion  of  Eversole  at  al,.  (26) 
that  the  diffusion  coefficient  nay  be  a function  of 
the  square  root  of  molar  concentration,  the 
linear  regressions  v:ore  repeated,  this  time  with 
X ® (Cq)^  and  Y as  before.  Now,  the  P-test  of 
the  data  shovxs  correlation  betv;cen  X and  Y at  the 
0,95  confidence  level.  The  upper  and  lov;er  bounds 
on  the  slope  of  the  regression  line  do  not  cross 
the  axis,  thus  indicating  a real  concentration 
effect.  The  best  line  through  the  data  is: 

Y ® -7#191X  + 6.944,  vrith  the  upper  and  lo’wer 


bounds  on  the  slope  at:  -1.839  and  respectively. 

(The  square-root  relationship  could  not  be  demonstrated 
with  the  35.5°  data;  this  is  probably  due  to  greater 
scatter  of  data  and  a smaller  number  of  points  to 
evaluate).'  As  shovm  in  Figure  8,  the  diffusion 
coefficient  of  lovz-conccntrations  of  CuSOi|,  in 
H2O  is  a v/cak  function  of  concentration.  This 
means  that,  for  lev;  concentrations,  it  may  be 
justifiable  to  neglect  concentration  effects  v;hen 
approximating  the  diffusion  coefficient.  As  stated 
earlier,  in  the  discussion  of  the  diffusion  equation, 
it  is  assumed  in  Pick's  First  Law  of  Diffusion, 
as  treated  in  this  study  (see  section  III),  that  the 
diffusion  coefficient  is  constant  v/ith  respect  to 
time,  distance,  and  concentration.  These  experiments 
indicate  that  that  assumption  is  justified  in  this 
range  of  concentrations. 

The  assumption  of  constant  diffusivity  with 
concentration  allov:s  the  calculation  of  a mean 
diffusion  coefficient  of  all  data  points  at  each 
temperature.  Of  the  38  experiments  at  25.0°C, 
in  the  concentration  range  0.005 ^ C© ^ 0.10,  a 
mean  diffusion  coefficient  of  5.^3  + 0.59  x 10"^  cm^/sec 
was  calculated  (955^  confidence  intfcrval  given) . 


In  the  same  concentration  range,  the  23  data  points 
taken  at  35.5'^C  lead  to  a mean  diffusion  coefficient 
of  6,675  i 0,97  X 10"^  cn^/sec.  The  double-tailed 
t-test  of  significance  chovrs,  even  at  the  0,93  level 
of  confidence,'  that  the  tv;o  values  are  different, 
suggesting  that  there  is  indeed  a temperature 
dependence  of  the  diffusion  coefficient,  as  predicted 
by  the  Stokes-Einstein  equation, 

Pigiu'e  9 is  a plot  of  -In  {^)  against  temperature 
for  the  tv;o  mean  values  discussed  above.  From  this 
plot,  an  activation  energy  of  diffusion  may  be 
estimated,  according  to  the  Boltzmann  Equation: 

0=  A exp  (-Ed/HT)  (6,3) 

where,  diffusion  coefficient,  cra^/sec 
A s=  frequency  factor,  cm2 /sec 

® activation  energy  for  diffusion,  cal/mole 
R «=  gas  constant  (1,987  cal/mole-OK) 

T = temperature,  °K, 

The  slope  of  the  bold  line  in  Figure  9 corresponds  to 
an  activation. energy  of  0,967  Kcal/mole,  The  dashed 
lines  represent  the  upper  and  lovior  bounds  of  the 
95^  confidence  limits  on  each  value,  corresponding 
to  values  of  Ed  of,  at  v;orst,  zero  at  minimum  and 
2,096  Kcal/mole  maximum.  These  values  are  in  the 
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same  order  of  macnitude  as  values  published  for 
similar  binary  systems  (20,24,29), 

An  interesting  aspect  of  the  temperature  effect 
vfas  noted  during  the  initial  experimentation  v/ith 
0,0511  6USO4  in  sulfuric  acid.  For  some  unknovm  reason, 
the  data  taken  at  that  concentration  led  to  the 
calculation  of  a diffusion  coefficient  much  higher 
than  that  of  the  other  CuS04  concentrations  in  the 
same  series.  Investigation  into  the  possible 
cause  of  this  obvious  error  turned  up  the  fact  that 
during  the  v;eek  in  which  the  0,05M  series  was  rim, 
the  air  conditioning  had  been  turned  off  every  night 
when  the  building  v;as  empty,  (this  v;as  during  the 
hottest  part  of  the  summer)  and  that  the  reservoir 
temperature  had  climbed  almost  to  30°C  overnight. 

When  the  air  conditioning  system  vrais  returned  to  I 

normal  24-*nour  operation,  the  0,05M  series  vias 
repeated,  and  the  initial  series  (with  no  temperature 
control)  vras  discarded,' 

B,  Stirred-Reservoir  Experiments 

In  order  to  evaluate  the  hydrodjuamic  effects 
caused  by  stirring  of  the  reservoir  material,  a series 
of  experiments  v/as  conducted  in  which  an  aqueous 
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solution  of  0.05  molar  CuSO^  diffused  into  flovjinn 
water  at  26,0  ± 1,0  ®C,  The  experiments  v/erc  divided 
into  tv;o  groups,  as  shovm  in  Table  5.5,  vfith  one 
series  run  at  a Reynold's  Number  (Wjjg)  of  2,03, 
and  the  other  at  Njj^  = 7.55.  In  each  series, 
a group  of  diffusion  experiments  vrere  run  at 
various  diffusion  times,  from  tv;o  to  tv;cnty-four 
hours  in  length.  The  tabulated  data  represent 
the  mean  and  standard  deviation  of  each  time  and 
experiment. 

For  each  group  of  values  at  a given  flow  rate, 
a linear  regression  v;as  run  so  as  to  determine  the 
effect  of  time  on  the  diffusion  coefficient.  That 
is,  it  v/as  desired  to  Imov;  if  diffusion  time  affected 
. the  apparent  diffusion  coefficient,  and  if  so,  at 
what  length  of  time  the  apparent  diffusion  coefficient 
approached  the  value  obtained  in  no-stirring  experiments. 
Using  the  variables  X = time  and  Y = diffusivity  x 10^, 
a linear  regression  of  tVie  data  at  = 2 led  to 
the  equation:  Y = -0,06X  + 6,831,  The  mean  value 
of  ^ for  the  whole  series  (and  95,^  confidence  limits) 

Is  6,1^6  i 0,761  X 10"^  cia2/sec.  The  upper  and  lower 
bounds  on  the  slope  of  the  regression  line,  respectively, 
are  +0,025  and  -0,l45,'  indicating  that  the  diffusivity 


docs  not  chance  with  diffusion  time.  Similarly, 
in  the  series  at  “Re  “ 7 .55,  shov?  no  effect  of  time 
on  diffusivity,  as  the  resrcssion  line  takes  the 
form  y = -0.077X  + 6.679.  Tliis  line  is  almost 
identical  to  the  regression  line  for  = 2.03. 

A double-tailed  t-test  of  the  two  sets  of  data  at 
two  Reynold’s  Numbers  indicates  no  difference 
between  the  tvro  groups  at  the  P = 0.05  level.' 

In  the  no-stirring  experiments  v;ith  0.05M  CuSO^, 
in  ivater  at  26®C,  the  mean  value  of  V obtained 
was  4.923  - .968  X 10**^  cm^/sec.  The  t-test 
indicates  that  the  mean  value  of  all  stirred- 
reservoir  experiments,  regardless  of  time  or 
flow  rate  (=  6. 01 5 x 10- <3)  is  not  different 
from  the  no-stir  value.  Interestingly,  hov^ever, 
if  only  the  short-time  (2  & 3 hour)  experiments 
are  considered,  at  both  flow  rates,  and  compared 
against  the  no-stir  experiment,  there  a 
difference  between  the  two  at  the  95^  confidence 
level.  On  the  other  hand,  comparison  of  the 
two  longest-time  experiments  (18  & 24  houi's) 
with  the  no-stir  experiment  does  not  shov;  any 
significant  difference. 

Thus  it  may  be  inferred  from  the  data  that 


for  experiments  in  v/hich  the  reservoir  is  stirred, 
longer  diffusion  times  approach  the  values  obtained 
in  no-stirring  experiments  better  than  the  short 
time  experiments.  If  the  so-called  "Ai-effect"  is 
leant,  then  its  effect  j-iay  become  negligible 
if  a long  enough  diffusion  time  is  allov/ed.  Hov/ever, 
as  long  as  a reliable  technique  exists  for  using 
the  no-stirring  model,  as  has  been  shovm  in  part  A 
of  this  section,  the  stirring  model  seems  impractical 
for  engineering  use. 

C.  Error  Analysis 

Possible  source^  of  error  in  the  capillary- 
reservoir  technique  include:  temperature  variation, 
concentration  measurement  inaccuracy.  Inexact 
filling  of  the  reservoir. and/or  capillary, 
measurements  of  time  and  capillary  length  and, 
for  the  stlrred-reservoir  system,  the  -O'^.-effect 
and  fluctuation  in  Reynold's  Number. 

Temi>erature  fluctuation  leads  to  an  error 
which  may  bo  approximated  by  Boltzmann's 
Equation  (equation  6.3),  in  vihich  ^ is  related  to 
l/T,  where  T is  in  degrees  Kelvin; 

A exp  (-%/RT) 
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Using  values  of  T = 2980K,  H = 1.9C7  cal/gnunole-OK, 
Ed  = 1000  cal/ginnole,  and  a maximum  temperature 
deviation  of  0,5®C  (this  vjas  the  worst  deviation 
recorded  in  any  no-stirring  experiment),  the 
maximum  error  in  '.'i  is  given  by: 

= 0'(E(i/o't)/RT2  = 0.011  S'*  (6.4) 


Thus  the  maximum  error  due  to  temperature  fluctuation 
is  about  1,1  in  the  no-stirring  series  of  experiments. 
For  the  stirred-reservoir  experiments,  a variation 
in  temperature  of  + 1®G  v;as  noted,  so  that  the 
'maximum  error  in  that  series  v.'ould  bo  about  2,2^, 


The  equations  in  section  III  xvhich  relate  iL. 
to  the  final-to-initial  concentration  ratio, 
all  have  the  general  form:  X = f(<i*  )* 

Thus  the  error  in  //  due  to  the  analytical  system 
may  be  expressed  as: 


<5Zf 
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1 ^9 

(6.5) 


Fi*om  the  definition  of  of  7^  , the  first  tv;o  terms  of 
equation  6,5  are  simplified  and  the  equation  becomes: 

C6-0 
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The  last  term  of  6,6  will  depend  on  the  equation 
for  ^ = f (X)  used.  The  Jt/t  terra  may  be  neglected, 
08  diffusion  time  is  carefully  controlled  - in  a 
typical  24-hour  run,  a one  minute  error  (the 


worst  that  miGht  be  expected)  only  contributes 
0#067J^  to  the  total  error.  The  2oL/L  term  is 
probably  not  neGlici^le,  and  is  estimated  by 
Luk  (20)  for  the  no-stirring  capillaries  to 
be 9 at  most,  around  3^^.  For  the  newly-made 
1-cm,  capillaries  used  in  the  stirred-reservoir 
experiments,  the  maximum  error  (based  on  Luk's 
formulation,  with  L = 1,0  rather  than  L = 2.54), 
is  8^.  The  last  term  of  equation  6,6  is 
determined  by  the  choice  of  equation  to  be 
used  in  the  concentration  analysis.  For  the 
no-stirred  experiments,  equation  3,11  v;as 
used: 

* 1 - (^)®’^^  (3.11) 

which  may  be  rearranged  as: 

(1  - (6,7) 

and  the  last  terra  of  equation  6,6  becomes: 

1-1^  . (1  - ^ 

Using  lowest  values  for  short-capillary  no-stirring 

experiments  (t  * 24  hrs.,'t‘=  0,235,  ^ = 0,481)  and 

allowing  for  the  precision  of  the  A,A.S,  at  ± IjS, 

S<fi  e 0,05,  and  the  error  term  comes  to  0,177, 

VO 


For  the  ctirred-reservoir  system,  Luk  (20) 
suggets  a "typical"  error  of  about  S^,  The 
inaxiraura  error  in  short  time  stirred  experiments 
(based  on  equation  3.9) » calculated  in  the  same  manner 
as  above,"  is  16,'^.  For  long-time  stirred-reservoir 
•experiments,  the  maximum  error  due  to  0 comes 
out  to  8^, 

Thus,'  summing  all  maximal  errors  for 
each  series,  the  follovring  are  given: 

No-Stirring:  max,  error  = 21,8^ 

Stirring  (trO.l):  max.  error  = 18,2^ 

^ r\  r\£\  . — — oli  nef  * 

In  all  experiments  discussed  earlier  in  this  section, 
only  the  no-stir  ZnO/KOH  experiments  showed  scatter 
beyond  the  maximum  error  range.  A possible  explanation 
of  why  this  happened  may  be  in  the  handling  of  the 
40/»  KOH  solution.  Because  of  the  danger  of  handling 
concentrated  base,  special  clothing  and  gloves  were 
worn  vdien  running  these’  experiments.  It  is  quite 
possible,  and  in  the  viev;  of  the  author,  even 
likely,  that  the  capillary  v/ashing  and  dilution 
procedure  was  not  as  carefully  folio-wed  as  in  the 
other  experiments,  ;ifhlch  required  no  special 
clothing  or  gloves.  Perhaps,  in  a situation  such 
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as  that,  v/hich  requires  special  handling,  many  data 
points  should  be  taken,  so  as  to  minimize  the  effect 
of  random  error  in  analytical  technique. 
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VII.  Conclur.lons 


!•  The  no-stirring  model  of  the  capillary-reservoir 
technique  gives  values  of  diffusion  coefficients 
in  line  with  those  predicted  in  the  literature  for 
the  three  systems  studied.  It  appears  to  be  a useful 
method  of  evaluating  diffusion  coefficients  of  these 
and  similar  binary  electrolytes, 

2.  The  stirred- reservoir  model  of  the  capillary 
technique  has  inherent  error  due  to  a hydro- 
dynamic  (A^)  effect.  With  long  diffusion  time 
■ the  effect  is  minimized. 

3*  Given  the  tv;o  conclusions  above,  the  no-stirring 
model  would  seem  to  be  a better  choice  for  a 
standard  laboratory  procedure  to  obtain  the 
integral  dlffusivities  of  binary  electrolytic 
8<^utlons,  at  a given  temperature  and  pressure. 

k.  Analysis  of  the  variance  of  data  in  the  CuS0|^/H2S0j^ 
system  showed  that  the  diffusion  coefficient  may  be 
a weak  function  (possibly  proportional  to  the  square 


root  of  Co  , as  suggested  by  others),  tut  that  this 
dependence  may  be  neglected  at  low  concentrations 


■ of  CuSOi,  such  as  were  studied  in  these  experiments. 
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The  Intenral  diffusivity  of  aqueous  IM  KCl  diffusing 
In  water  at  25^0  was  found  to  be  1.843  ± 0,113  x 10-5  cn^/sec. 
This  value  is  in  good  agreement  v;ith  published  data 
for  the  same  system. 


6,  The  integral  diffusivity  of  O.OIH  ZnO  in  40^  KOH  ] 

diffusing  in  40J»  KOK  solution  was  measured  as 

7.22  ± 1,90  X 10-6  cm2 /sec  using  the  no-stir  i 

capillary-reservoir  technique.  Because  of  the 

large  scatter  and  small  number  of  data  points, 

no  significant  difference  could  be  found  betv;een 

this  value  and  that  of  a similar  system  using 

23^  KOH  solution,  although  the  Stokes-Einstein  ■ 

equation  predicts  a difference  due  to  viscosity 

effects.  For  this  reason  it  is  suggested  that 

future  users  of  this  technique  allov;  a sufficient 

I 

number  of  data  points  to  minimize  scatter, 

7.  The  integral  diffusivity  of  lov;  concentration 
CuS0|^  in  half-molar  sulfuric  acid,  diffusing 

In  0,5H  vras  found  to  have  a significant  ^ 

temperature  dependence.  At  25.0®C,  the  value 
is  5.43  i 0.59  X 10-^  cm2/sec.  At  35.5°C, 
across  the  same  range  of  concentrations,  the 
Integral  diffusivity  is  6,675  x 10“6  cm2 /sec, 

TOio  range  of  concentrations  studied  was  0, 005-0, lOM.  i 

i 
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8,  The  CuS0/^/H20  system  was  studied  using  both  the 
stirring  and  no-stirring  models  of  the  capillary- 
reservoir  method*  The  no-stir  experimental  results 
gave  a value  of  the  integral  diffusion  coefficient 
as  4,923  i 0,968  X lO"^  cn^/sec,  in  good  agreement 
with  published  data  for  the  same  system.  The 

i • flow  channel  experiments,  v;hsn  averaged  over  a 

I 

i range  of  diffusion  times  from  tv7o  to  twenty-four 

! hours,  did  not  shovi  significant  difference  from 

' the  no-stir  values  (v;ithin  95/^  confidence  limits). 

However,  when  only  diffusion  times  of  three  hours 
or  less  are  considered,  the  average  diffusivity  is 
significantly  higher  than  that  found  in  the  no-stir 
experiments.  It  is  hypothesized  that  this  is  due  to 
a hydrodynamic  effect  which  increases  the  .apparent 
diffusion  coefficient  due  to  convection  of  diffusant 
out  of  the  capillary.  This  effect  is  minimized 
with  long  diffusion  time. 

9. '  The  activation  energy  for  diffusion  in  the  low 

concentration  CuSOji^/HgSO^  system  studied  was  found 
to  be  on  the  order  of  1 Kcal/gm-mole, 


1 


VIII.  Recomnandntions 


1, '  The  effect  of  viscosity  and  pH  on  the  diffusion 

coefficient  iray  be  studied  in  the  ZnO/KOH  or 
CuS02^/H2S0|^  system  as  follovjs:  Choose  a 
concentration  of  primary  electrolyte  (e.g.,  ,05M  ZnO) 
and  a range  of  supporting  electrolyte  concentrations 
0-4M  KOH,  in  l-molar  increments),  end  run 
several  no-stir  diffusion  experiments  at  each 
ooncentraticn  of  supporting  electrolyte.  Then, 
repeat  the  series  of  experiments  but  viithout 
the  supprting  acid  or  base  electrolyte.  Instead, 
use  a nantral-TtH  v^fscous  «r»T*nRfi!  « •nannre* 

of  viscosities  comparable  to  the  initial  series 
of  experiments.  Comparison  of  the  tv/o  series 
of  experiments  v/ill  indicate  how  the  apparent 
diffusivity  is  affected  by  acid/base  concentration, 
and  by  viscosity. 

2, -  It  would  be  interesting  to  study  the  feasibility 

of  using  the  capillary-reservoir  technique  on 
binary  systems  other  than  electrolytes,  such 
as,  for  example,  toluene/benxene  mixtures, 
biological  fluids,  etc. 
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Symbol  Heanlnt 


Inside  diameter  of  capillary,  cra,‘ 
time,  hr,*  or  sec,'  (as  noted) 
width  of  flow  channel,  cm, 

t 

length  vector:  distance  from  bottom  of  capillary 

horizontal  distance  vector,  cm,‘ 

vertical  distance  vector,  cm, 

cross-sectional  area  of  flow  channel,  cra2 

concentration  of  diffusant,  moles/liter 

initial  "C"  at  start  of  experiment 

average  cone,  of  diffusant  in  capillary  after 
allotted  diffusion  time, 

depth  of  water  in  flov;  channel,  cm, 

capillary  length,  cm, 

effective  diffusion  length,  cm, 

Reynold’s  Number 

volumetric  flow  rate,  ml, /sec, 

temperature,  degrees  Celcius 

velocity,  cm,/sec, 

centerline  velocity,  cm, /see,' 

surface  velocity  above  centerline,  cm, /sec, 

diffusion  coefficient,  cm  /sec. 


non-diffusing  region  of  capillary  caused  by  flow 


w > 


7k  . 


A ' 

density  at  specified  temp,’  and  pressure,  gm/sra-^  I 

I 

viscosity  at  " ■ ■ « ^ poise  | 

dimensionless  time  (see  p.’  20)  I 

dimensionless  concentration  ratio  (see  p,  20)  \ 

I 

dimensionless  length  (see  p.  20)  I 

rotation  rate  ^ 
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APPENDIX  B:  Calculation  of  Or>tinm  Capillary  Leri.rth 

There  are  interdependent  factors  to  be  considered 
when  choosing  the  length  of  the  capillary  diffusion  coll; 
namely,  diffusion  time  and  depth.  At  the  beginning  of 
the  diffusion  time  pcrios,  the  concentration  of  diffusant 
is  uniform  (C  = C^)  throughout  the  cell.  After  immersion 
in  the  stirred  reservoir,  the  interfacial  concentration 
goes  to  zero  and  a concentration  gradient  of  diffusant 
develops  along  the  length  of  the  capillary.  As  the 
diffusion  progresses,  the  depth  of  the  gradient  increases, 
until  eventually  it  touches  the  closed  end  (bottom)  of 
tile  capillary. 

As  explained  earlier,  the  diffusion  coefficient 
is  computed  on  the  basis  of  measurement  of  diffusant 
remaining  in  the  capillary  after  a measured  time  period. 
Therefore,  error  in  diffusant  concentration  measurement 
results  in  error  in  calculated  diffusivity.  It  has  been 
shown  (27)  that  the  error  in  calculated  diffusivity 
associated  v/ith  a given  concentration  measurement 
error  is  highly  dependent  on  the  amount  of  diffusant 
remaining  in  the  capillary.  For  example,  if  90/i  of 
the  initial  amount  of  diffusant  present  in  the  capillary 
remains  there  after  diffusion,  a Ifi  error  in  concentration 
measurement  corresponds  to  1 17,7^  error  in  the  resultant 


calculated  diffusivity.  In  a diffusion  cell  in  v;hich 
50^  of  the  initially-present  diffusant  remains, 
the  same  error  in  concentration  measurement 
corresponds  to  a mere  2%  error  in  calculated  diffusivity,' 
In  order  to  reduce  this  to  about  1%,  over  70/»  of  the 
initially-present  diffusant  would  have  to  be  removed,' 

The  stirring  experiments  were  to  be  run  for  a 
broad  range  of  diffusion  times,  so  as  to  determine  the 
magnitude  of  the  iW-region  as  a function  of  time.  It 
was  considered  desirable  to  have  diffusion  times  no  longer 
than  a day,  y\:t  allovf  no  worse  than  a 10^  error  in 
calculated  diffusivity,  (The  arjalytical  method 
of  analyzing  capillary  contents  has  an  inherent  1% 
uncertainty.)  In  such  a case,  a tv;o-hour  experiment, 
a one-centimeter  long  capillary  would  just  meet  this 
requirement, 

A capillary  shorter  than  one  cm,’  would  not  be 
practical,  as  a too-small  capillary  volume  leads  to  a 
great  error  in  diffusant  removal.  That  is,  the  larger 
the  capillary  volume,  the  less  significant  a spilled 
drop  appears  in  the  analysis.  Thus,  the  optimum 
practical  capillary  length  for  stirred-reservoir 
experiments,  1 mm,  capillary  diameter,  with  0,05  Molar 
CuS04  in  water,  would  be  one  centimeter. 
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APPENDIX  C:  Sample  CalculatlOT?3 


1.  Calibration  of  tho  A.A.S. 


Calibration  solutions  of  CuSO^  in  distilled  H2O 
were  Trade  ,up,  as  shovm  bolov;,  and  were  run  on  the 
A.A.S. , and  readings  (on  an  arbitrary  :;ero-to-100 
scale)  v.'ero  recorded: 


Cone,  of  Cu^  ::  IQ^M 

3 

4 

5 

' 6 


A.A.S.  Rcadin'T 

36.0 

46.0 
■ 56.5 

67.0 


Letting  "C"  be  concentration  x 1o5h,  and  "AA" 

the  reading  on  tlio  A.A.S.,  a least-squares  reduction 

line  through  the  four  points  above  gives  the  equation: 

AA  * 10.35" 0.  +4.80 


This  calibration  equation  v;as  used  to  obtain  concentration 
values  of  ion  in  sol-ution  from  readings  on  the  A.A.S. 


ZJ  ExnT.r)Ie  of  no-stirri>v  rr::''ol 

One  no-stirring  experiment  which  v.’as  run  at 
the  same  tine  as  the  above  calibration,  was  v;ith 
0.05M  CuSO;^  in  0.51-  H2S0^^^  at  J!5°C.  There  were  9 
data  points  taken,  but  only  one  is  presented  here 


to  illustrate  the  calculations.  One  such  point 
was  for  a 3«^9cn,  long  capillary,  initially 
filled  with  0.05fi  CuSOji^/acid  solution,  which 
had  diffused  for  24  hours  in  0.5I'  acid,  Tlie 
A,A,S,  reading  for  this  capillary,  vxhich  had  been 
v;ashed,  and  its  contents  diluted  to  50  ml,,  was 
AA  = 64,0,  From  the  calibration  equation  on  the 
previous  page,  the  concentration  of  copper  in 
the  diluted  sample  was  C = 5.7199  x 10“5molar, 

To  obtain  a value  of  (^,  the  firial-to- initial 
concentration  ratio,  divide  C by  capillary  volume 
(0,0701  ml),  multiply  by  dilution  factor  (50  ml), 
and  divide  by  initial  concentration  (0,05M): 


4)  *=  0,816 

As  this  is  a short-time,  no-stir  expci'iment, 
use  equation  3.11  to  calculate  a dimensionless  time: 


e, . 1 


(3.11) 


Rearrangement  of  the  equation  so  as  to  solve 
for  T gives  a value  of  T=  0,035.  From  the 
definition  of  T , we  obtain: 

fY  - _ I.  ;>  / 


I 


4,919  X 10“^  cm^/scc 


‘ i 


3.‘  The  stirred-reservoir  experiments 

I ■ . 

I The  same  calibration  method  vras  used  in  those 

It 

I experiments  as  for  the  no-stirring  series.  Likewise, 

the  calculation  of  concentration  ratio,  (|),  was  made 
in  the  same  fashion.  The  choice  of  equation  to 
. ’ use  to  calculate  T was  based  on  the  value  of  ^ 

j calculated.  If  ({)  >0.72  (corresponding  to  T = 0.06) , 

< the  short-time  approximation,  equation  3*9  vras 

used.'  If  9^^0.633  (corresponding  to  T = 0.10),  the 
long-time  approximation,  equation  3.8,  iras  used. 

For  values  of  ({)  betv^een  the  limits  described  above,  ! 

' a five-terra  expansion  of  the  series  solution  (eq.  3,6)  ; 

I 

was  solved  by  Newton's  Method.'  A computer  program  j 

was  written  i-xhich  calculated  ({)»  chose  the  appropriate 
equation,  calculated  t , and  from  that,  computed  . 

Example:  For  one  capillary,  length  1 cm.,  ! 

diffusion  time  = 24  hours,  with  0.05H  CuSOjij,  in  , 

i 

water  at  26oc,  <{)  was  calculated  as  as  0.246.  This 
falls  in  the  category  of  the  long-time  approximation, 

I 

so  the  equation  used  ^ras  3.8: 

I 

9 “ ^ • exp  (-TT2r/4)  (3.11) 


jk 


. Rearrangement  of  the  equation,  solving  for  X.  , gives 
a value  of  ^ = 0.484,  or  5.602  x 10”^  cm^/scc. 


no  ^ 


k.'  Analysis  of  Varlar.co  of  Dp.ta 

The  ZnO/KOH  data  discussed  in  section  VI-A-2 
will  be  used  as  an  example  of  the  statistical 
analysis  of  data.  The  basic  data  are  found  in 
Table  5.2,  on  page  45. ‘ A prepared  program  for 
use  in  a Hewlett-Faclcard  9100A  programmable  desk 
calculator,  program  09100-70801,  calculated  the 
mean  and  standard  deviation  of  each  of  two 
input  samples,  "X",  and  "Y",  and  computed  the 
t-statistic  for  correlation  of  the  two  samples. 


For  example,  the  program  vjas  used  to  compute 


O 


Table  5.2.  The  95%  confidence  limits  of  the  data  * 
were  calculated  by  multiplying  the  computed  standard 
deviation  by  the  tabulated  t-statistic  for  7 degrees 
of  freedom  and  P=0,05  (30),  V/here  N = number  of 
data  points  ( =8),  t = 2,365,  and  the  computed 
standard  deviation  is  2,09  x 10-^,  the  calculated 
95%  confidence  interval  is: 


0*9^  * s t 2,42  X 10-6  . 

The  99^  confidence  interval  \«is  computed  in  the  same 
manner,  except  that  the  t-statistic  used  v;as  for 
P •*  0,01  rather  than  0,05,  When  one  data  point  was 


} 


I 


>1 


found  outside  the  interval,'  it  was  discarded 
and  a nev/  mean  ^ireis  calculated,* 

To  compare  tvio  samples  of  data,  the  prosrara 
described  above  vias  used  to  calculate  a t-statistic 
for  comparison  ;>rith  tabulated  ■values.  For  example, 
in  compai’ing  the  data  of  Table  5.2  in  terms  of 
capillary  length,  the  data  vzas  divided  into  two 
groups:  sample  "X"  included  points  taken  v;ith 

2.  *5^  cm  capillaries  and  sample  “Y"  ^iias  the 

3. ^9  cm  capillary  sample.  Tlie  program  computed  a 
t-statistic  of  0.6l  for  the  two  samples,  far  less 
than  the  tabulated  value  of  2,57  (for  5 degrees 

. of  freedom) , thus  indicating  no  difference 
betv/een  the  two  samples,  and,  no  proof  of 
any  effect  of  capillary  length  on  the  diffusion 
coefficient  for  these  data  points. 
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